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1 Introduction 

In their paper "Should All Crystals Be bcc? Landau Theory of Solidification 
and Crystal Nucleation" [1] , S. Alexander and J. McTague review the Landau 
theory of solidification [2] which predicts that solids in three space dimensions 
form a bcc(body centered cubic) structure and pointed out that indeed un- 
der certain circumstances many solids near their melting point exhibit this 
behavior. Landau's theory of phase-transitions implies that an octahedron 
structure in momentum space gives a global extremum to the free-energy, 
and thus is the preferred structure. This is true in three spatial dimensions. 
In two dimensions either triangular or honeycombed lattices are predicted to 
form depending on the sign of the appropriate order parameter. Solidification 
occurs when an inhomogeneous configuration becomes energetically favored 
compared to a homogeneous one. The condensing density waves allow vari- 
ous patterns of spontaneously broken translational and rotational invariance. 
For a review of somewhat more complex spontaneous breaking of space time 
symmetries see i.e. [5]. 

In string theory higher dimensional systems are natural, in fact ten dimen- 
sional supersymmetrical systems have a special role as being stable. Assuming 
the microscopic theory allows the formation of solids in higher dimensions it 
is interesting to investigate if general structures can emerge generically. In 
tachyonic string backgrounds closed and open tachyon fields are also forced 
to have a non-zero wave number, condensing this field leads to structures sim- 
ilar to that of solidification, however not necessarily in three dimensions [1]. 
This motivates generalizing the study of Alexander and McTague to higher 
dimensions, which is the subject of this paper. Following their work we as- 
sume the structure of regular polytops and find that the result remains true 
in higher fiat non compact dimension. 

In section 2 we review the work of S. Alexander and J. McTague and in 
section 3 we generalize it to more than three space dimensions. 
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2 Review 



In Landau's theory of phase transition, one expands the free-energy, 

$ = $2 + ^3 + - (1) 

Where $„ is of order n in the density field p. 

The presence of a $3 term is essential for the solidification to occur. A $4 
term provides stabilization. 

The second-order term, in momentum-space, is: 

^2 = j d?qid?q2A {qt, ql) S {qt + g|) = J d^qA (q) p-^p_-^ (2) 

With pg- being the Fourier component of the density function, i.e. 
p(r) — 10 ^3/2 / d^qPqC''^^"^ . The form of the function A{q) depends on the 
microscopic theory. The delta-function ensures translational-invariance and 
that A{q) is a function of a single q. Moreover isotropy implies that A {q) 
depends actually only on the magnitude of q. Translational invariance will 
be spontaneously broken for those class of functions A (q) for which the mag- 
nitude Q of the vector q does not vanish at the extremum. In such a ground 
state the configuration pg gets support only on the momentum space sphere 
of this radius. The second order term is then proportional to / dQ^p^p_-'Q, 
which equals to / d^rp^{r)^ henceforth denoted by p^.The magnitude of p^ is 
fixed by the fourth order term. 

In momentum-space the third-order term is: 

^3 = / d\d^q2d%B {qt, -^,ql,T)S{qt + -^ + ql) Pf^p-^p-^ (3) 

Here T stand for all of the thermodynamical quantities on which B depends. 
Since the magnitude of the gt's is already fixed, one can re-write this term 
as: 

$3 = Bq (t) j dn^dn^dn^s (ot + + oJ) p^^p^/^^ (4) 

One sees that $3 gets a contribution only from configurations in which the 
— yi 

Qi s form equilateral triangles. Thus rotational invariance is also sponta- 
neously broken. Since the position-space density (p(r)) is real, each configu- 
ration has to satisfy the equation P ^ — P'^^ 1 "^^^ich gives us 
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as possible density configurations. For such a configuration, with 2n values 
of q for which pq ^ 0, (if q is such a value, then necessarily -q is also such 
a value), in order to maximize $3, at a given p^, all p-^f's should be of the 

2 

same magnitude. One obtains that n 
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= p^, so the term Pj^ Pj^ Pj^ is 

VI V2 V3 



proportional to n ' . One can use this result to determine the configuration 
of the crystal. The momenta Qi for which p^f 7^ form some polyhedron 

in momentum-space. This polyhedron has exactly 2n edges. The $3 term 
gets a contribution from each face of that polyhedron which is an equilateral 

-if-] ■ Nt, where Ne is the number of 
edges of the polyhedron, and Nf is the number of faces that are equilateral 
triangles. The regular, convex, three-dimensional polyhedron that maximizes 
the $3 term, is the octahedron [Ij. Although two tetrahedra give the same 
result, this is essentially the same configuration, i.e. the momentum vectors 
that generate one generate the other. 



3 Higher-Dimensional Space 

In higher dimensions, the same ideas apply. One considers the regular, convex 
polytopes in higher dimensions. Postponing for a moment the discussion of 
the case of a four-dimensional space, which is a bit more complicated, we 
discuss five or more dimensions. It is known [5] that in five dimensions or 
more there are only three regular polytopes, all convex: The n-simplex, the 
n-hypercube, and its dual the n-cross polytope, which is the n-dimensional 
analog of the octahedron. The free energy corresponding to each polytope 
configuration depends on the number of elements in its one-skeleton (edges), 
and on the number of elements of its two-skeleton (two-dimesional faces) that 
are equilateral triangles. The hypercube does not contribute in this case, since 
its two-skeleton does not include any triangles, equilateral or otherwise, so 
the value of $3 for it is zero. For the n-simplex the number of elements 
in its one-skeleton is given by ("^^^ |S] and the number of elements in its 

two-skeleton is given by ("3^) [5]. All the elements of the two-skeleton are 
equilateral triangles, so we should consider all of them in our calculation. 
One should notice, though, that a single n-simplex does not contain any 
pair of opposite edges, so two n-simplexes have to be considered, which gives 
Ne = 2 ■ ("^^) elements in the one-skeleton and Nt = 2 ■ elements in 
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the two-skeleton. The $3 term is then proportional to 

Turning our attention now to the cross-polytope, we have 2^ ■ elements 

in the one-skeleton [5], and 2^ ■ elements in the two-skeleton [5], all of 
which are equilateral triangles. This will give a value of 
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3-y/n(n — 1) 



Dividing these two results, we get that the ratio between the value of $3 for 
the cross-polytope and the value of $3 for the simplex is 

n — 1 \n — 1/ 

which is greater than one for n greater than three. So we see that for five 
dimensions, or more, the cross-polytope, which is the n-dimensional analogue 
of the octahedron, is the preferred momentum-space configuration in these 
dimensions as well. In the special case of four dimensions one has six convex 
regular polytopes [5]. A direct check shows that here also the cross-polytope, 
called a 16-cell in four dimensions, is preferred. Thus the result holds for any 
dimension greater than two. Since ([7]) is a rising function of n for n > 2, a 
higher dimensional polytope will always be preferable to a lower dimensional 
one, ruling out liquid-crystals as a possible structure for a high-dimensional 
crystal, for the case studied here of a scalar order parameters. Tensorial order 
parameters allow more structure already at n = 3. To determine the structure 
of the resulting lattice in position-space, notice that the vertices of the n cross 
polytope are of the form Cj, 1 < i < n, with ±1 in the i'th coordinate, and 
in all other. This means that the edges of this polytope are n-vectors, with 
±1 in the i'th and j'th place, i 7^ j, and in all others, i.e. a vector of the 
form 

0,0,..., ±^ ,0,0,..., ±^ ,0,...,0 (9) 

y i'th place j'th place J 

Thus, the lattice in momentum-space is generated by n linearly-independent 
such vectors. The reciprocal lattice, in position-space, is the set of vectors 
where scalar products with the momentum lattice vectors are integers. It 
is generated by vectors with ±| in all their entries. This is exactly a b.c.c. 
lattice. 
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4 Conclusions 



We have shown that the result that near melting a bcc lattice is preferred, is 
true in any number of flat non compact dimensions. 

Several issues require further study. One is how to generalize this analysis 
to compact dimensions of various topologies. The global structure of the man- 
ifold will further constrain the allowed configurations requiring in some cases 
to deal with commensurability aspects. Another is the case of configurations 
that form polytopes that are not convex or not regular, as may occur on the 
sphere. It is possible that such a configuration, which would be preferable to 
the cross-polytope, exists. Such a configuration will not form a crystal, but 
it is possible, if it has sufficient symmetry, that it will form a quasi-crystal. 
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